RM L52K14a 


NACA 

RESEARCH MEMORANDUM 

AN ANALYSIS OF THE LATERAL STABILITY 

OF THE DOUGLAS D-558-H AIRPLANE EQUIPPED WITH A YAW 

DAMPER, WITH SPECIAL REFERENCE TO THE EFFECT OF 

YAW-DAMPER RATE -GYRO SPIN-AXIS ORIENTATION 

By Ordway B. Gates, Jr., Albert A. Schy, 
and C. H. Woodling 

Langley Aeronautical Laboratory 
Langley Field, Va. 


NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 

WASHINGTON 

March 16, 1953 
Declassified May 8, 1957 


NACA RM L52Kl4a 


NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


RESEARCH MEMORANDUM 


AN ANALYSIS OF THE LATERAL STABILITY 
OF THE DOUGLAS D-558-II AIRPLANE EQUIPPED WITH A YAW 
DAMPER, WITH SPECIAL REFERENCE TO THE EFFECT OF 
YAW-DAMPER RATE-GYRO SPIN-AXIS ORIENTATION 

By Ordway B. Gates, Jr.,, Albert A. Schy, 
and C. H. Woodling 

SUMMARY 

A theoretical investigation has been made to determine the effect of 
a. yaw damper on the lateral stability of the Dougla.s D-558-II airplane 
for various anticipated flight conditions. Since the airplane angle of 
attack varies considerably with changes in flight condition, the incli- 
nation of the rate-gyro spin axis is very important in determining the 
stabilizing effect of the yaw damper. Possible advantages of orienting 
the gyro spin axis in the direction of the airplane Y-axis a.re considered. 
A combination of yaw-damper gearing ratio and yaw-damper axes inclination 
is selected so that the airplane has satisfactory Dutch roll stability 
for the flight conditions considered. These flight conditions were 
cruising at a. Mach number of 1.6 at 5G,000 feet and 70,000 feet in 
lg flight, at the same altitudes and Ma.ch number in 2g flight, and 
landing at two different lift coefficients. 

A particular auxiliary control surface is recommended for use with 
the d amp er, and theoretical estimates are presented of the aerodynamic 
characteristics of the proposed surface. The results of the investiga- 
tion indicated that a combination of gearing ratio and yaw-damper inclina- 
tion can be selected so that the lateral stability characteristics of the 
Douglas D- 558 -II should be satisfactory in all the flight' conditions 
considered. 


INTRODUCTION 


Recent flight tests of the Douglas D-558-II research airplane have 
indicated that the lateral or Dutch. roll oscillation of this airplane is 
very poorly damped. At present, the undesirable stability character- 
istics of the airplane in the cruising or clean condition at supersonic 
Mach numbers are of primary interest since these characteristics must be 
improved before the airplane . can reach the speeds of which it is poten- 
tially capable. Also, some difficulty has been experienced in the 
landing condition. A yaw damper has been designed and constructed by 
the Flight Research Division of the Langley Aeronautical Laboratory and, 
when installed in the airplane, it should markedly improve the poor 
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Dutch roll stability of the Douglas D-558-II. Briefly, this system will 
cause an auxiliary surface to be deflected proportional to the airplane 
rate of yaw in such a manner as to damp the lateral motions of the air- 
plane. A previous analysis of the effects of a yaw damper on the sta- 
bility of the Douglas D-558-II has been reported in reference 1, but the 
present investigation makes use of stability derivatives obtained from 
recent wind-tunnel tests at both subsonic and supersonic speeds and, in 
addition, takes into consideration the actual characteristics of the yaw 
damper . 

The purpose of this paper is to determine by theoretical analysis 
the effects of the proposed yaw damper on the lateral stability of the 
Douglas D-558-II for several representative flight conditions of the air- 
plane and to recommend, on the basis of these effects, a gearing ratio 
(control deflection per unit rate of yaw) and an inclination of the rate 
gyro Z-axis to the airplane Z-axis such that the poor damping of the air- 
plane lateral oscillation will be satisfactorily improved. The term 
"satisfactory damping," as used in this paper, implies that the lateral 
oscillation damps to one-half amplitude in one cycle or less. This is 
an arbitrary value which was chosen for convenience in the calculations 
and is somewhat more rigid than the Air Force criterion (ref. 2) for 
satisfactory damping over the range of frequencies of the Douglas D- 558 -II 
In addition, attention is given to the design and location of the control 
surface which is to be actuated by the yaw damper, since the surface 
characteristics will have a direct bearing on the selection of a suitable 
gearing ratio for the yaw damper. 

The results of this analysis are presented in the form of stability 
boundaries and airplane motions for the various flight conditions. The 
motions were calculated on a Reeves Electronic Analog Computer. 


SYMBOIE AND COEFFICIENTS 

0 angle of roll, radians 

\|r angle of yaw, radians 

p angle of sideslip, v/V, radians 

r, ijr yawing angular velocity, dty/dt, radians/sec 

p,0 rollixg angular velocity, d0/dt, radians/sec 

Q pitching angular velocity, d0/dt, radians/sec 

v sideslip velocity along lateral axis, ft/sec, radians/sec 
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airspeed, ft/sec 
Mach number 

pitching moment of gyro about gyro Y-axis 
mass density of air, slugs/cu ft 
dynamic pressure, | pV 2 , lb/ sq ft 

wing span, ft 

wing area, sq ft 

weight of airplane, lb 

mass of airplane, w/g, slugs 

acceleration due to gravity, ft/sec 2 

relative -density factor, m/ pSb 

angle between longitudinal body axis and principal 
axis, positive when body axis is above principal 
axis at the nose, deg 

inclination of principal longitudinal axis of airplane 
with respect to flight path, positive when principal 
axis is above flight path at nose, a - e, deg 

angle of flight path to horizontal axis, positive in a 
climb, deg 

airplane moment of inertia about principal X-axis, 

p O 

I x = mk x , slug-ft^ 

airplane moment of inertia about principal Z-axis, 

Iz = mkz Q 2 , slug-ft 2 

gyro moment of inertia about gyro Y- or gimbal axis 


gyro moment of inertia about gyro X- or spin axis 
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radius of gyration in roll' about principal longitudinal 
axis, ft 

radius of gyration in yaw about principal vertical 
axis, ft 

nondimens ional radius of gyration in roll about prin- 
cipal longitudinal axis, jb 

nondimensional radius of gyration in yaw about prin- 
cipal vertical axis, jb 

nondimensional radius of gyration in roll about longi- 
tudinal stability axis, \j*c 2 cos^t] + ^sin^t] 


nondimensional radius of gyration in yaw about vertical 
stability axis, 


\f. 


o p p p 

COS T) + Ky Sin T] 
o 1 A o 


nondimensional product-of-inertia parameter, 

3in T] COS T) 


(^o 2 ' K X 0 2 ) Sil 


trim lift coefficient. 


W cos 7 


qS 


rolling-moment coefficient, 
yawing-moment coefficient, 
lateral-force coefficient, 


Rolling moment 
qSb 

Yawing moment 
qSb 

Lateral force 
qS 
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,n r n rb 
d 2V 


'"p“ 


l V 

2V 


' y p a 

2V 


^ rb 
d 2V 


time, sec 

differential operator, d/dt 

period of oscillation, sec 

time for amplitude of oscillation or aperiodic mode to 
change by factor of 2j positive value indicates 
damping, negative value indicates divergence 

cycles to damp to one-half amplitude, 

real part of complex root of characteristic stability 
equation 



cn 


angular frequency, radians/sec 
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yaw-damper natural frequency, radians/sec 
yaw-damper damping ratio 


deflection of the auxiliary control surface, radians 

deflection of pump servo arm (see fig. l) 

dC n dC 2 

control effectiveness parameters, - — ; r — 

bS A ’ bb A ■ 

angle of attack of longitudinal body axis to flight 
path, deg (see fig. 2) 

rate -gyro spin axis 

1 

rate -gyro gimbal axis 

rate-gyro axis normal to spin and gimbal axes 

inclination of rate-gyro Zg-axis to the airplane 
Z-axis, deg (see fig. 2) 

angle between rate -gyro Z^-axis and flight-path Z-axis 
when airplane is undisturbed, deg (see fig. 2) 

angular displacement of rate-gyro spin axis about gyro 
Y- or gimbal axis due to airplane angular velocity 
about rate-gyro Z-axis, deg (see fig. 2) 

inclination of rate-gyro Z-axis to flight-path Z-axis; 
(k = k q + Sq when gyro spin axis is oriented in the 

direction of airplane X-axis; k = k q when oriented 
in the direction of airplane Y-axis ), deg 

airplane acceleration about gyro Y- or gimbal axis, , 
radians/sec^ 

yawing velocity about airplane Z-axis, 

\J/ S cos a 0 + 0 g sin a Q , radians/sec 
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yawing velocity about an axis perpendicular to the gyro 
spin axis, ^g cos k + 0g A sin k, radians/sec 


^pi 


M= 


gyro precessional moment due to airplane rate of yaw 
about gyro Zg-axis 

gyro restoring moment due to gyro deflection, dMy/dSg 


% 


M, 


gyro damping moment due to rate-of-gyro deflection, 

SMy/S&G 

moment about gyro gimbal axis due to Coulomb friction 






rotational velocity of gyro about gyro spin axis, 
radians/ sec 


gearing ratio of yaw damper, 


b A /i G 


K 


gearing ratio, 


& A/ & G 


servo time constant 


l 

b 


nondimensional distance from center of gravity of air- 
plane to center of pressure of auxiliary control 
surface 


h 

b 


nondimensional distance from airplane longitudinal body 
axis to center of pressure of auxiliary control 
surface 


Subscripts: 

G gyro axis 

S stability axis 

B airplane body axis 


A 


auxiliary control surface 
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DESCRIPTION OF THE PROPOSED YAW DAMPER 
AND AUXILIARY SURFACE 
Yaw Damper 


The proposed yaw damper consists of a spring -re strained gyro which . 
is deflected proportional to the airplane rate of yaw about the gyro 
Z-axis and a hydraulic .servoac tuat or which deflects an auxiliary control 
surface proportional to the rate-gyro spin-axis deflection. A sketch of 
the system is presented in figure 1 and the system of axes employed in 
the analysis is shown in figure 2. Briefly, when the airplane experi- 
ences a rate of yaw about the gyro Z G -axis, a precessional moment Mp r is 

applied to the gyro which is equal to the product of the gyro moment of 
inertia about its spin axis, the angular velocity of the gyro, and the 
airplane rate of yaw about the gyro Z-axis; that is, 

Mpr = %%% (1) 


The differential equation which describes the gyro spin-axis deflection 
subsequent to an applied moment is 


Iy G D 2 - M 5 G D - M S g ) 5 G = ^ " si S n 1)6 X M f (2) 


where I v is the moment of inertia of the gyro about its Y-axis, Mg 
y G c 

is a measure of the gyro viscous damping, Mg is the spring constant 

G 


of the restraining springs, ' Mf is the moment due to Coulomb friction 

(the notation "sign D5" is used to point out that this moment always 
opposes the gyro rate of deflection), % is any applied moment about 

the gyro Y-axis, which for the problem being considered is that moment 
defined by equation (l), and b is the airplane angular acceleration 
about the gyro Y-axis; hence, 




!X 


" I Y g ® " si § n DS X M f 


( 3 ) 


The moment due to Coulomb friction Mf has been measured and found to 
be approximately 0-33 inch-pound. The minimum value of ijr G , based on 
this value of M f , below which the system will be insensitive is about 
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1/4° per second. Since this value is small, M f is assumed to be zero 

throughout the analysis to follow. The angular velocity of the airplane 
about the gyro Z-axis is related to the airplane yawing and rolling 
velocities about stability axes by the following expression: 

= ^s cos K + 0s sin K (*0 


where k is the angle between the stability Z-axis and the gyro Z-axis. 

If the gyro spin axis is oriented in the direction of the airplane 
X-axis, the angle k will vary, since the gyro is deflected about its 
gimbal axis (gyro Y-axis) due to airplane yawing about the gyro Z-axis. 

If, however, the gyro spin axis is oriented in the direction of the air- 
plane Y-axis, the angle k is essentially independent of the gyro deflec- 
tion) hence, in the discussion to follow both orientations are considered. 

Gyro spin axis oriented in direction of airplane X-axis (fig. 2(a)) . - 
For this orientation, k may be defined as follows: 

k = k . + 8 P 

o Cr 

where is the angle between the gyro Z-axis and the stability Z-axis 

when = 0, and Sq is the gyro deflection about its Y- or gimbal 

axis due to airplane yawing about the gyro Z-axis. The angle k q varies 

from one flight condition to another since it is dependent upon the air- 
plane trim angle of attack. The angularity between the gyro Z-axis and 
the airplane-body Z-axis (subsequently referred to as $>) is set on the 
ground aiiu does not vary as the flight condition is changed; thus, k q 

is defined as 


k „ = a. 


0 


Then, k becomes 




and, for the assumption that the cosine and sine of a small angle are 
approximately equal to unity and to the angle in radians, respectively. 


= t s + ( a o - °)^s + Ms 


( 5 ) 
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Thus, when the gyro spin axis is oriented in the direction of the air- 
plane X-axis, is a nonlinear function of and jZS g due to the 

term 5 q 0 s • 

Gyro spin axis oriented in direction of airplane Y-axis (fig. 2(b)) 
If the spin- axis of the rate gyro is oriented so that it is coincident 
with the airplane Y-axis when = 0 and if the gyro Z-axis is tilted 

at an angle 0 to the airplane Z-axis, then the angular velocity of the 
airplane about the gyro Z-axis is given as 


* G = ^ [q] cos 5 G + 0 sin (6) 

^G ® 

where \Jr run is the yawing velocity of the airplane about an axis 

LJ &g=° 

which is coincident with the gyro Z-axis when Sq = 0 and 9 is the 

airplane pitching velocity. The expression for r_,-i is 

L G J 5q=0 

i[ G ] &G=0 = cos K o + K sin K o (7) 

Substituting this expression into equation (6) and making the usual 
simplifications for small angles gives the equation 

i G = + (a Q - $)0 S + es G (8) 

which is similar to equation (5) except that, for this orientation, the 
nonlinearity is due to 9 rather than 0 S . Since, in the derivation of 
the equations of motion, lateral and longitudinal motions are considered 
to be independent and 9 generally is much less than 0 B in lateral 
maneuvers, the assumption is made that = 0. Thus, throughout the 

remainder of this paper the rate -gyro spin axis is considered to be 
oriented in the direction of the airplane Y-axis. If a future investi- 
gation should indicate that the nonlinear term 8 g 0 s has a negligible 

effect, the results presented- in this paper are applicable to either 
gyro orientation. 
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Substitution of equation (8) into equation (3) yields the result 
( for 05 q = 0 and a = $ G ) 

(ly G D . 2 - M 6 G D - M 6 g ) 5 G = Z X G %^s + ( a o - *)h\ - I Y g ^G 

(9) 

where '0 G = '0 S - (a Q - <t») V Since I Yq £ I Xq , = 8l6 radians per 

second, and |0 G | a: lo|\j/ G |, the term ly^G should contribute negligibly 
to the results of this analysis and is omitted. 

Relation of gyro-spin-axis deflection to auxiliary- surface deflec- 
tion.- The deflection of the auxiliary control surface as a function of 
the gyro-spin-axis deflection & G and the hydraulic servo time con- 
stant t is given by 

(1 + t D)6 a = K5 g (10) 

The constant K, in terms of the various lengths indicated in figure 1, 
is 


K = 


ni z 
m + n wy 


(id 


The time constant of this servo has been found to have a negligible 
effect on the results of this analysis and, hence, the assumption is 
made that S A = K6 q. For this assumption, equation (9) may be expressed 

as 


(d 2 + 2|cd 0 D + cd 0 2 )& a = K 0 o) 0 2 % + (a Q 


M P 


where av 


d G . 


is the damping ratio of the rate 


and K q , which is referred to as the static gearing 

KIv CDy 
X G X G 



( 12 ) 


gyro 

ratio 
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For the values of n, i, 


m. 


v. 


y> 


in table I, the gearing ratio K 0 is, 


V 

as a function of 


°v 


and Me 


given 


z, K Q = 2.66 z. 


The' range over which z may be varied is from 0 to 3.2 inches; hence, 

K 0 is variable from 0 to approximately 8.5* The maximum deflection of 

the pump servo arm 5 pg (fig. l) has been measured and found to be ±22.5° 


This deflection is, in terms of the gyro-arm deflection, = — & p . 

’ P B mh u 

For the values of n, i, m, and h given, the pump servo arm is 
deflected to its maximum value when & G = ±3.24°. In order to prevent 

damage to the linkages, stops should be placed on the gyro to limit the 
gyro deflection to ±3-24°. From equation (3), it can be seen that & Q 

steady state will be equal to 3.24° for \jr G = 0.125 radian per second. 
Therefore, this yaw damper can be expected to operate linearly only if 
ty(j is less than this value. It should be pointed out, however, that, 
if the maximum deflection of the control is assumed to be ±20°, the con- 
trol would be deflected to its maximum value for ty G less than 
0.125 radian per second if K Q is greater than 2.8. Thus, for higher 
values of K Q , the maximum value of ^ G for which the yaw damper could 

be expected to operate linearly would be less than 0.125 radian per 
second and, hence, for K 0 >2.8 the physical limits of the control 

surface would be more restrictive than the limits on the gyro-spin-axis 
deflection. 


Auxiliary Surface 

The auxiliary surface, when deflected, introduces a yawing moment 
about the airplane Z-axis and a rolling moment about the airplane X-axis; 
that is. 


( C “)b C »6 a 5a 

( c 4 ■ V a 


The effect of the sideforce introduced by the surface is assumed to 
be negligible except in the calculation of the yawing and rolling moments 
produced by the deflection of the surface. The yawing and rolling 
moments about the stability or flight-path system of axes are related to 
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the moments about body axes as follows: 


( c 4 = ( c 4 cos a ° ■ ( C *) B sln a ° 

( c j) g ■ ( c l) B cos a o + ( c n) B sln a o 


or, for the assumption of small angles of attack, 

Co, • ('■*.), • f 




> (13) 


The auxiliary surface proposed for use with the yaw damper is shown 


in figure 3* The surface is to be mounted 13 feet forward of and 2^- feet 

6 

below the airplane center of gravity on the undersurface of the fuselage, 
slightly offset to the left of center (looking forward from center of 
gravity). This offset is necessary in order that the surface not inter- 
fere with the lowering and retraction of the nose wheel. The surface is 
to be hinged at approximately the midpoint of the root chord, and counter- 
clockwise rotation of the surface about its hinge, as seen looking down, 
is taken as positive rotation; thus, positive deflection of the surface 
results in a negative sideforce or a negative yawing moment about the 
airplane Z-axis. This convention is selected in order that K Q be posi- 
tive for a yaw damper which introduces a yawing moment opposing the air- 
plane yawing velocity. The values of C 


n? 


and C 


h 


of this surface 


related to airplane body axes were estimated from the results presented 
in references 3 and ^ and unpublished data to be -0.01 and 0.0022 per 
radian, respectively. These values were found to be essentially the 
same for all flight conditions considered. The values of C n and . 

b a 5 a 

about stability axes can be determined from equation ( 13 ) for each flight 
condition by substituting the airplane angle of attack and are presented 
in table II. 
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EQUATIONS OF MOTION 


The equations which describe the lateral motion of the yaw-damper- 
equipped airplane about stability or flight-path axes are as follows: 


Yawing, 


L„Kz 2 S§ D 2 


- 1 s ? d )+ + ( 2M -b K xz $ d2 - 1 % i 




Rolling, 


b 2 J2. 


i 2 Ub K XZ ^2 D 


l Clr | D U + *i b K x 2 ^D 2 -|u z 


- % Cl | DW 
V c 2 P v / 


> (14) 


C 2„P = C Z R 8 A + C Z 


Sideslip, 


Control, 


(a^b | D )* - c l(* + (2% I D - ^p) 13 = 0 


^D 2 + 2 |oj 0 D + % 2 j& A - K 0 o> 0 2 D n|r + (a Q - ®)D0j 


= 0 


The stability derivatives C Y , C Y , and C Y have been neglected 

5 a P r 

in equations (14); also, the equations are for level flight. These- equa- 
tions are assumed to represent the case for which the rate-gyro spin axis 
is oriented in the direction of the airplane Y-axis, and the gyro Z-axis 
is at an angle <J> to the airplane Z-axis in the airplane X,Z plane. 

The. characteristic equation of the airplane yaw-damper system is of 
the form 


AX 


'54^2 
BX 2 + CX -+ DX° + EX 


+ FX + G = 0 


(15) 
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where the coefficients A, B , C, and so forth for any flight condition 
of the airplane and given values for the yaw-damper parameters £ and co Q 

can be expressed as functions of the yaw-damper gearing ratio K Q and 
the gyro Z-axis inclination 0. 


ANALYSIS AND DISCUSSION 


Effect of linear yaw damper .- Calculations were made for four 
cruising conditions at a Mach number equal to 1.6 and two landing condi-. 
tions of the Douglas D-558-H in an effort to select a yaw-damper gearing 
ratio K 0 and an inclination of the gyro Z-axis to the fuselage Z-axis 
such that the airplane has satisfactory Dutch roll stability for all the 
flight conditions considered. The mass and aerodynamic characteristics 
for each of the investigated flight conditions are presented in table II 
and were estimated from the results of references 5 and 6 and unpublished 
wind-tunnel tests at supersonic Mach numbers. For each of flight condi- 
tions 1 to 5 (see table II) a boundary was obtained in the K Q , 0-plane 

which defined the combinations of K Q and $ for which the Dutch roll 

oscillation would damp to one -half amplitude in one period of the oscil- 
lation or less, and the results are presented in figure k. These bounda- 
ries were calculated by methods similar to those presented .in refer- 
ences 7 and 8. The cross-hatched region shown in the figure defines the 
combinations of Kq and $ for which the Dutch roll oscillation of all 
the flight conditions will damp to one-half amplitude in less than one 
period, but the values of Kq within this region are so large as to 
make the region impractical. Satisfactory oscillatory stability can be 
obtained for flight conditions 1 to 4, however, for much smaller values 
of Kq. Calculations were made on the Reeves Electronic Analog Computer 
for cases 1 to k in which Kq and 0 were varied systematically and 
for K Q = 2.5 and 0=2° all four cases were very satisfactory. This 

combination of gearing and inclination was selected on the basis of 
satisfactory stability and low gearing ratio and because slight varia- 
tions from these values produce small variation in the system stability. 
For case 5 and 0 = 2°, a gearing ratio of approximately 6.5 would be 
necessary for the oscillation to damp to one-half amplitude in one cycle. 
However, for K Q = 2.5 and 0 = 2°, the oscillation of case 5 damps to 
one-half amplitude in about 1.5 cycles which satisfies the existing Air 
Force criterion (ref. 2) for period -damping relationship of the lateral 
oscillation. Thus, if this is considered satisfactory for the landing 
case, no problem arises. The yaw-damper gearing ratio can be adjusted 

in the cockpit from 0 to about 8.5; hence, if V\/2 “1*5 is not 
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acceptable for the landing case, the value of Kq can be increased to 
about 6. 5 by the pilot and the stability for this case should then be 
satisfactory. 

The results of the calculations made on the Analog Computer are 
presented in figure 5 for all six flight conditions for Kq = 2.5 and 
$ = 2°. Also, for case 5> motions are shown for K Q = 6-5 and 0=2°. 

For comparison, the motions for the respective flight conditions without 
yaw damper are also shown, and in every case except case 5 the increase 
in stability is very apparent. Flight condition 6 is presented, since 
this condition has been used in flight tests as a means of determining 
the characteristics of the D-558-II in the landing configuration. The 
lateral motions for this case without yaw damper are seen to be less 
stable than those obtained for the true landing case (case 5) , hut the 
yaw damper is far more effective in the simulated landing configuration 
than in the true landing case. The increased effectiveness of the yaw 
damper for the landing condition at = O.36 (case 6) as compared 

with that obtained for C L = 1.05 niay he attributed, to a large degree, 
to the higher flight speed for C L = O.36. This conclusion is based on 

the fact that this yaw damper essentially increases the damping- in- yaw 
stability derivative which generally has a very appreciable stabi- 

lizing effect' on the stability of the Dutch roll oscillation. The incre 
ment to this derivative, if the lags in the yaw damper are negligible, 
can be shown to be 


2C 




2Ko c % , 


V 

b 


Thus, this increment to is seen to vary directly with the air- 

speed V and, for the landing configuration at = O.36 at 

12,000 feet, would be approximately twice that obtained for the landing 
case at C L = 1.05 at sea level or 12,000 feet. It should be pointed 

out that in these calculations the control deflections were limited to 
±20°, since this is assumed to be the physical limit on the auxiliary- 
surface deflection. As was pointed out previously, the gyro-arm deflec- 
tion is limited to 3. 2k° which occurs for = +0.125 radian per second 

Hence, the limits on the control deflection for K Q = 2.5 should have 

been approximately ±18°, but this slight difference should have a negli- 
gible effect on the results presented. 

Additional calculations were made to determine the variation of 
period and damping in the vicinity of K Q = 2.5 and 0=2° for all 
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six flight conditions and the results are presented in table III. These 
results are in agreement with those obtained from the previously dis- 
cussed Analog studies, and for every case except the landing condition 
at sea level, the Dutch roll oscillation should damp to one-half ampli- 
tude in close to, or less than one cycle throughout the range of K Q 

and <t> considered. The stability characteristics, presented in table III, 
of the airplane in the landing condition at sea level, although not as 
desirable as those of the other flight conditions considered, might be 
considered satisfactory. The primary conclusion which should be reached 
from the results shown in table III is that values of K Q = 2.5 and 

$ = 2°, selected as the optimum combination of control gearing and gyro 
Z-axis inclination for the investigated flight conditions as a whole, 
are not critical values since the airplane stability for combinations 
of K 0 and $ in the vicinity of the chosen values is not appreciably 

different from that calculated for K Q = 2.5 and 0=2°. 

Effect of variations in the principal moments of inertia .- The 

values of K Y 2 and Kn 2 from which were calculated K x 2 , Kg 2 , 

A o o 

and for the various flight conditions were obtained from data pre- 

sented in reference 5. More recent data on the airplane principal 
moments of inertia which were provided by the NACA High-Speed Flight 
Research Station at Edwards Air Force Base, Calif, (referred to herein 
as HSFRS) , differed considerably from those of reference 5 and for pur- 
poses of comparison both sets of values are presented in figure 6 plotted 
against wing loading w/S. Throughout the calculations, the assumption 
is made that the principal longitudinal axis is located 3-7° below the 
fuselage axis and is invariant with wing loading. Calculations were made 
for each of flight conditions 1 to 6 in which the moments of inertia 
afforded by HSFRS were utilized. The results of these calculations indi- 
cated no appreciable effect of the assumed variations in the inertial 
moments and are not presented, but it may be of interest to note that the 
largest effects were obtained for the landing configurations. 


SUGGESTIONS FOR FUTURE RESEARCH 


A more detailed analysis of the effects of the various nonlinearities 
in the yaw-damper system, such as Coulomb friction, operation of the 
control surface at full deflection, and pitching of the gyro spin axis 
following airplane rotational velocities, should be made. Also, the 
effect of these nonlinearities on the airplane response characteristics 
following various inputs should be investigated. 
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CONCLUDING REMARKS 


An analysis made to determine the effect of a yav damper, designed 
and constructed by the Flight Research Division, on the Dutch roll sta- 
bility of the Douglas D-558-II research airplane indicates that this yaw 
damper, when used in conjunction with the auxiliary control surface 
described in this paper, should be capable of improving markedly the 
inherently poor stability of the airplane. Calculations made for several 
representative flight conditions indicate that satisfactory .stability 
should be obtained for the high-speed flight conditions for a control 
gearing of 2.5 when the gyro gimbal axis is inclined 2° nose down rela- 
tive to the airplane X-axis, but the stability is marginal for the 
landing configuration at sea level. If better stability is desired for 
the landing case a gearing ratio of 6.5 should be used for this gyro- 
axis inclination. 


A preliminary phase of the analysis showed that the rate-gyro spin- 
axis deflection varied nonlinearly with airplane rate of roll about the 
stability X-axis if the gyro spin axis was oriented along the direction 
of the airplane X-axis but varied linearly with rate of roll if the spin 
axis was oriented along the airplane Y-axis . For this reason, the gyro 
spin axis was assumed to be oriented along the airplane Y-axis for all 
the calculations presented in this paper. 


If the control gearing is to be variable in the cockpit from 0 
to 8.5., stops should be placed on the gyro to limit its deflection to 


1 ° 


approximately ±3— 


so as to prevent damage to the linkages. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va. 
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TAB IE I 

YAW-DAMPER PARAMETERS 


z X n > slug-ft 2 . . 
I v , slug-ft 2 . . 
Mg , ft-lbs/radian 


Mg , ft-lbs/radian/sec 


Yaw-damper lever-arm lengths, in. : 

z 

n 

i 

m 

w 

x 

y 

h 

ojy , radians/sec 

X G 


I 

cd 0 , radians/sec . 


O.OO765 

0.0090 

-13.8 

-O.386 


variable 
. . 5.0 

. . 5-0 

. . 1.8 
. 0.625 

. . 3.0 

. . 1 
. 2 

. . 8l6 
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TABLE II 

STABILITY DERIVATIVES AND MASS CHARACTERISTICS USED IN 
CALCULATIONS OF LATERAL STABILITY OF D- 558 -II 




Flight condition and configuration 



1 

Clean 

II 

Clean 

111 

Clean 

IV 

Clean 

V 

Landing 

VI 

Landing 

l/w 

1 

0.2 

1 

0.2 

1 

1 

C L 

0.473 

0.095 

0.184 

0.037 

1.05 

0.36 



707 

707 

275 

275 

33 

47.5 

e, deg .... 

3.7 

3.7 

3-7 

3.7 

3.7 

3.7 

a Q , deg ... . 

6.6 

0.28 

1.8 

- 0.7 

- 9.7 

- 1-7 

t), deg .... 

2.9 

-3.42 

- 1.9 

- 4.4 

6 

-5.4 

K x 2 

0 . 0157 H 

0.015833 

0.015534 

O.O16116 

0.016729 

0.016772 

K z 2 

0 . 13669 

0.13657 

0 . 13687 

0 . 13628 

0.13567 

0.13563 

>Scz 

0.0061141 

-0.0072056 

- 0. 0040093 

-0.0092559 

0 . 012579 

-0.011336 

Cll P 

C, 

0.087 

0.087 

0.087 

0.087 

0.23 

0.29 

- 0.035 

-0.067 

- 0.057 

-0.074 

- 0.20 

-O.15 

P 

c Y 

Y P 

N 

-0.726 

-0.726 

-0.726 

-0.726 

- 1.15 

-1.43 

- 0.50 

-0.54 

- 0.51 

-0.56 

-0.63 

-0.56 

% 

0.076 

0.143 

0.122 

0.159 

0.13 

0.17 



0 

0 

0 

Q 

n 

0 

% 

% 

-0.045 

- 0.0084 

- 0.017 

-o.oo 4 

-O.17 

-0.06 

-0.25 

- 0.25 

- 0.25 

-0.25 

- 0.33 

-0.31 

° y p 

0 

0 

0 

0 

0 

0 

tan 7 

0 

0 

0 

b 

0 

0 

V, ft/sec . . . 

1553 

1553 

1553 

1553 

225 

460 

S, ft 2 .... 

175 

175 

175 

175 

175 

175 

b, ft 

25 

25 

25 

25 

1 25 

■ 25 

W, lbs .... 

14,000 

14,000 

14,000 

14,000 

11,000 

11,000 

Alt., ft . . . 

70,000 

70,000 

50,000 

50,000 

0 

12,000 

Mach number . . 

1.6 

1.6 

1.6 

1.6 

0.20 

0.43 

w/s 

80 

80 

80 

80 

63 

| • 63 

Ns ■■■■ 

-0.01 

-0.01 

-0.01 

-0.01 

- 0.01 

- 0.01 

Ns ■ • • • 

0.001 

0.0022 

0.0019 

0.0023 

0.0005 

0.0025 
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TABIE III 


VARIATION OF PERIOD AND DAMPING WITH $ AND Kq FOR 
THE VARIOUS FLIGHT CONDITIONS CONSIDERED 


K o 

Aperiodic Modes 

Oscillatory Modes 

T l/2 

T l/2 

T l/2 

P 


(a) Case 1 


0 

89.4 

1 . 8 o 

13.5 

4.30 

3.14 

2.5 

12.1 

1.49 

3-40 

4.43 

• 77 




• 033 

.197 

.17 

2.5 

11.5 

1.56 

3.16 

4.44 

• 71 




.033 

• 197 

• 17 

2.5 

10.9 

1.65 

2.95 

4.45 

.66 




.033 

• 197 

• 17 

2.5 

10.3 

1.74 

2-77 

4.46 

.62 




.033 

.197 

.17 

2.0 

13.0 

1.68 

3.54 

4.43 

.80 




.033 

• 197 

• 17 

2.5 

10.9 

1.65 

2.95. 

4.45 

.66 



1 . 6 l 

.033 

• 197 

.17 

3.0 

9.4 

2.53 

4.48 

• 57 




.033 

.197 

.17 


(b) Case 2 


0 

41.90 

1.08 

-7.30 

5-51 

- 1.30 

2.5 

4.74 

1.16 

7-39 

5.80 

1.27 




• 033 

• 197 

• 17 

2.5 

4.04 

1.33 

5.56 

5 .'91 

• 94 




• 033 

• 197 

.17 

2.5 

3.17 

1.64 

4.38 

6.01 

• 73 




.033 

.197 

• 17 

2.5 


— 

3.54 

6.11 

• 58 




2.40 

61.42 

.04 




.033 

.197 

' .17 

2.0 

4.76 

1.38 

6.6 3 

5.86 

1.13 




.033 

.197 

• 17 

2.5 

3.17 

1.64 

4.38 

6.01 

• 73 




.033 

• 197 

• 17 

3.° 


— 

3.22 

6.20 

.52 




2.21 

51.96 

.043 




• 033 

• 197 

• 17 
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TABLE III 


VARIATION OF PERIOD AND DAMPING WITH <5 AND K Q FOR 
THE VARIOUS FLIGHT CONDITIONS CONSIDERED 
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TABLE III 


VARIATION OF IERIOD AND DAMPING WITH $ AND K Q FOR 
THE VARIOUS FLIGHT CONDITIONS CONSIDERED 






























26 


NACA RM L52Kl4a 


Z 


i 7 
0,6 

i 

l 



X 


B 


G 


0 


(a) Gyro spin axis oriented in direction of airplane X-axis. 



(b) Gyro spin axis oriented in direction of airplane Y-axis. 
Figure 2.- System of axes used in analysis. 


L-77031 


















Figure 6.- Variation of moments of inertia of the Douglas D-558-II with 

wing loading. 



